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High-harmonic generation from hydrogenlike muonic atoms exposed to ultraintense high-frequency 
laser fields is calculated. Systems of low nuclear charge number Z are considered where a nonrel- 
ativistic description applies. By comparing the radiative response for different isotopes we demon- 
strate characteristic signatures of the finite nuclear mass and size in the harmonic spectra. In 
particular, for Z > 1, an effective muon charge appears in the Schrodinger equation for the relative 
particle motion, which influences the position of the harmonic cutoff. Cutoff energies in the MeV 
domain can be achieved, offering prospects for the generation of uffrashort coherent 7-ray pulses. 
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I. INTRODUCTION 

One of the most successful and accurate methods to 
probe nuclear properties employs muonic atoms [l|, 
Due to the small Bohr radius of these exotic atoms, the 
muonic wave function has a large overlap vi^ith the bind- 
ing nucleus. Precision measurements of muonic tran- 
sitions to deeply bound states can therefore reveal nu- 
clear structure information such as finite size, deforma- 
tion, surface thickness, and polarization. The first x-ray 
spectroscopy of muonic atoms was performed in 1953 us- 
ing a 4- meter cyclotron Q. Today, large-scale facilities 
like TRIUMF (Vancouver, Canada) or PSI (Villingen, 
Switzerland) exist which are specialized in the efficient 
generation of muons and muonic atoms H. New de- 
velopments aim at the production of radioactive muonic 
isotopes for conducting spectroscopic studies on unstable 
nuclear species Q. Muons bound in atoms are also able 
to catalyze nuclear fission [6] and fusion 1'^ reactions. 

On a different front, the field of laser- nuclear physics 
is emerging [8|. While lasers have always represented 
important tools for nuclear spectroscopy Q, in recent 
years their role is changing qualitatively and growing 
because of the tremendous progress in high-power laser 
technology. The interaction of intense short laser pulses 
(/ 10^* — 10^°W/cm^) with matter can produce 
highly energetic electrons, protons and photons (e.g., via 
bremsstrahlung) . In pioneering experiments, this has led 
to the observation of laser-induced nuclear fission p^ . 
nuclear fusion |11|, and neutron production in nuclear 
reactions [l^- Advanced laser sources might also pave 
the way to nuclear quantum optics 11311 and coherent 7- 
spectroscopy using ultrashort pulses |14l - [l6| . 

In light of this, the combination of muonic atoms with 
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intense laser fields opens promising perspectives. Con- 
trary to the traditional spectroscopy of muon transi- 
tions between stationary bound states, the exposure of 
a muonic atom to a strong laser field renders the prob- 
lem explicitly time-dependent and the muon, thus, a dy- 
namic nuclear probe. In this setup, the muon is coher- 
ently driven across the nucleus which, for example, gives 
rise to the emission of radiation and, in general, allows for 
time-resolved studies on a femtosecond scale. The infor- 
mation on the nucleus gained by laser-assistance can in 
principle complement the knowledge obtained from the 
usual field-free spectroscopy of muonic atoms. 

Against this background, we have recently considered 
the process of high-harmonic generation (HHG) from 
strongly laser-driven muonic hydrogen and deuterium 
atoms p/7|] . The process of HHG represents a frequency 
up-conversion of the applied laser frequency due to a non- 
linear coupling of the atom with the driving external field 
(see [18l - [22| for recent reviews). It can be understood 
within a three-step model, where the bound lepton is 
liberated from the atom by tunneling ionization, propa- 
gates in the laser field, and finally recombines with the 
core, returning its kinetic energy upon photoemission. 
By way of a comparative study it was demonstrated that 
the harmonic response from muonic hy drog en isotopes is 
sensitive to the nuclear mass and size [13] • This shows 
that muonic atoms subject to strong laser fields can re- 
veal information on nuclear degrees of freedom. Muonic 
deuterium molecules in superintense laser fields represent 
another interesting example towards this combined effort, 
where field-induced modifications of muon-catalyzed fu- 
sion have been investigated ^] . Muonic hydrogen atoms 
have moreover been studied as systems which could allow 
for observation of the Unruh effect [23] . 

In this paper we extend our previous study on HHG 
p^l to hydrogenlike muonic atoms (ions) with nuclear 
charge number Z > 1. To this end, the time-dependent 
Schrodinger equation in the presence of the binding nu- 
cleus and a few-cycle laser pulse is solved on a numerical 
grid in one spatial dimension. From the resulting dipole 
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acceleration the HHG spectrum is extracted and com- 
pared for different isotopes. Characteristic effects aris- 
ing from the finite nuclear mass and size are revealed. 
Moreover, in the case Z > 1 the laser-driven particle 
dynamics becomes more complex because the center-of- 
mass of the atomic constituents does not stay at rest any 
longer. As a result, an effective muon charge appears in 
the Schrodinger equation for the relative motion, which 
affects the harmonic cutoff position. The cutoff energies 
achievable with muonic atoms in the nonrelativistic do- 
main of interaction are very large, reaching several MeV. 
This holds in principle prospects for the production of 
coherent 7-ray pulses of ultrashort duration. 

Since muonic atoms are tightly bound systems, laser 
fields of extraordinary field strength and photon en- 
ergy are required to infiuence the muon motion. In 
the ground state of muonic hydrogen, for example, the 
muon is bound by 2.5 keV and experiences a binding 
Coulomb field strength of 1.8 x 10^'' V/cm corresponding 
to the field intensity 4.2 x 10^^ W/cm^. A comparison 
of these numbers with the parameters of the most ad- 
vanced present-day and near-future laser sources is use- 
ful. In the range of optical and near-infrared frequencies 
{fkjj ^ 1 eV) , the highest intensity presently attainable is 
~ 10^^ W/cm^ [2^ and the next generation of high-power 
lasers aims at intensities of 10^'^W/cm^ and beyond 
[Hi. In the VUV frequency domain (fiw - 10 - 100 eV) 
a maximum intensity of ^ lO^^W/cm^ has been at- 
tained with a free-electron laser at the FLASH facil- 
ity (DESY, Germany) [13 . The Linac Coherent Light 
Source (SLAC, Stanford) has recently entered the fre- 
quency domain hoj ^ 1 keV [2^ . Near- future upgrades 
of such machines are planned to produce brillant x-ray 
beams (fiw ^ lOkeV) with peak intensities close to 
10^0 W/cm^. There are also efforts to generate ultra- 
short, high-frequency radiation {huj ^ 10 — 1000 eV) from 
plasma surface harmonics where considerably higher in- 
tensities might be reachable due to a high conversion effi- 
ciency [2^ . With these novel sources of intense coherent 
radiation it will become possible to influence the quan- 
tum dynamics of light muonic atoms with nuclear charge 
numbers Z < 10. Note that in principle the laser field 
strengths and frequencies experienced by the atoms can 
be enhanced further when a relativistic atomic beam is 
employed instead of a fixed target [13, dH . 

As to their lifetime, we point out that light muonic 
atoms and molecules may be regarded as quasistable sys- 
tems on the ultrashort time-scales of strong laser pulses 
(t ~ fs-ps), since their life time is determined by the free 
muon life time of 2.2 /is. For the field parameters assumed 
in this paper, the influence of the external laser field on 
the muon decay is immaterial as well [111- In deeply 
bound states of heavy atoms, the muon life time can be 
reduced due to absorption by the nucleus to ~ 10~®s 
which still exceeds typical laser pulse durations by or- 
ders of magnitude. 

We organize the paper as follows: Section II deals with 
the theoretical framework in which the separation into 



center-of-mass and relative coordinates of the two-body 
Schrodinger equation for a hydrogenlike muonic atom in a 
laser field is performed. We also give here a scaling trans- 
formation between ordinary and muonic atoms, as well as 
suitable model potentials which allow us to incorporate 
the nuclear mass and size. Section III has been reserved 
for the presentation of our numerical results and their 
discussion. While Sec. III. A compiles the harmonic cut- 
off energies available from different low-Z muonic atoms, 
Sects. III. B and III.C show a series of calculations de- 
voted to the impact of the nuclear mass and size on the 
HHG spectra. A comparison of the nuclear signatures 
predicted for muonic atoms with those to be expected in 
highly-charged electronic ions is undertaken in Sec. III.D. 
The conclusion of the paper is given in Sec. IV. 



II. THEORETICAL FRAMEWORK 

A. Separation of relative and center-of-mass 
motion 

We consider the nonrelativistic quantum dynamics 
of an initially bound muon in a few-cycle laser pulse 
described by the time-dependent Schrodinger equation 
(TDSE). For our laser parameters of interest, we may 
ignore the space dependence of the laser field (dipole 
approximation), treating it as a purely time-dependent 
electric field. Due to the large muon mass, the atomic 
nucleus cannot be considered as infinitely heavy. We 
therefore start from the two-particle TDSE written in 
the length gauge as: 
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where and to„ are the muonic and nuclear masses, x^ 
and the coordinate vectors for the muon and the 

nucleus, and — —iHd/dx^ and p„ = —iHd/dxn the 
corresponding momentum operators, respectively. Be- 
sides, the nuclear charge number is Z, the elementary 
charge unit e, the binding potential ^ (|x^ — x„|), and 
the laser electric field E(t) which oscillates with angular 
frequency a;. 

The application of the dipole approximation in Eq. ([T]) 
considerably simplifies the problem as it renders the 
Schrodinger equation ([T]) for the muon-nucleus two-body 
system separable into relative and center-of-mass motion. 
By introducing relative and center-of-mass coordinates 
x = x^ — x„ and X = (m^x^ -t- m„x„)/M, respectively, 
with the total mass M = rUf^ + m„, one finds that the 
evolution of the center-of-mass is described by 
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with P = —itld/d'X-. Equation ^ is the non-relativistic 
Volkov equation for a particle of charge {Z — l)e and mass 
M in the presence of a laser field. As a consequence, the 
center-of-mass motion does not emit higher harmonic fre- 
quencies and may therefore be ignored in the following. 
Note that in the special case Z — 1 (i.e., hydrogen iso- 
topes) the center-of-mass moves freely, while the laser 
field only couples to the relative coordinate 



The relative motion is governed by (see also |23l. l34j) 
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with the reduced mass — m^rrin/M , the relative mo- 
mentum p = —iHd/dx, and the effective charge 
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In the special case Z — 1, the effective charge reduces to 
Qe — e, whereas Qe > e holds for atomic numbers Z > 1. 

Formally, Eq. ^ is the Schrddinger equation for a sin- 
gle particle of charge —qe and mass rrir in the presence 
of a nucleus and a laser field. The accelerated motion 
of the relative coordinate in the combined external fields 
therefore gives rise to the emission of higher harmonics. 
We note that, in physical terms, the relative coordinate 
accounts for the fact that both the nucleus and the muon 
oscillate in the laser field with different amplitudes in op- 
posite directions (see also Fig. 2 in [iTj). 

We point out that the effective charge in Eq. Q has 
already been derived in Ref. [35|. It was also shown 
there that the two-body TDSE in a laser field separates 
straightforwardly in the velocity gauge. It is interesting 
to note that, in contrast, the separability in the length 
gauge is a more subtle issue because the operations of 
performing gauge transformations and dipole approxima- 
tions are not commutative (35l. [36|. As a consequence, 
cross terms appear in the exact version of the two-body 
TDSE in the length gauge which, strictly speaking, pre- 
vent the equation from being separable. The cross terms 
typically become important at the borderline to the rela- 
tivistic regime when the value of the relativistic field pa- 
rameter [see Eq. approaches unity [35|. For the non- 
relativistic laser parameters applied in the present study, 
however, these terms are very small and have therefore 
been neglected in Eq. (P). 



B. Scaling considerations 

The form of Eq. ([3]) is equivalent to the Schrddinger 
equation for an ordinary hydrogen atom (i.e., an elec- 
tron bound to an infinitely heavy proton) in a laser field. 
The relation can be made explicit by virtue of a general 
method known as scaling transformation. For the spe- 
cial case of a laser-driven atom, the scaling procedure is 
nicely explained in [sjj. Suppose, that we have a hydro- 
genlike muonic system of nuclear charge Z on the one 



side and an ordinary hydrogen atom on the other side. 
We introduce an electronic coordinate vector Xe and time 
te and relate them to the muonic coordinate x and time 
t of Eq. ^ according to 
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with the mass ratio p = rrie/mr and the electron mass 
me- When rewritten in the scaled space and time, Eq. ^ 
becomes 
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with = —ihd/dxe and the scaled laser frequency and 
field strength 
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This means that a muonic hydrogenlike atom in a laser 
field with parameters E and lu behaves like an ordinary 
hydrogen atom in a field with E,, and We given by Eq. ([7]) , 
provided that the binding potential V{x) arises from a 
pointlike nucleus. To give an example, the typical pa- 
rameters of an intense Ti:Sapphire laser Huje — 1.5 eV, 
Ee = 2.7x 10^ V/cm (10^^ W/cm^) translate to a muonic 
helium atom as Huj = 1.2 keV, = 9.1 x lO^^V/cm 
(1.1 X 10^^ W/cm^). This comparison demonstrates that, 
despite the huge laser intensities applied in our computa- 
tions, the laser-driven muon dynamics remains nonrela- 
tivistic due to the large muon mass. Moreover, since the 
laser intensities required for HHG from muonic atoms 
are very large already for hydrogen isotopes and steeply 
increase with the nuclear charge, we restrict our consid- 
eration to muonic atoms with Z < 10. An advantage 
of these systems as compared to heavier ones is that the 
relative differences in mass and size among isotopes are 
larger for \ow-Z atoms in the nuclear chart. 

We emphasize that the scaling procedure does not ac- 
count for nuclear properties like the finite nuclear size or 
the nuclear shape. Evidently, when the transition from, 
e.g., a muonic hydrogen atom to an ordinary hydrogen 
atom is performed, the proton radius is not to be length- 
scaled in accordance with Eq. ([5]) but remains fixed. As a 
consequence, for atomic systems where nuclear properties 
play a role, not all physical information can be obtained 
from the knowledge of the ordinary-atom case via scal- 
ing. In Sec. III.C below we show results which display 
the influence of the nuclear size on the process of HHG. 



C. ID approximation and model potentials 

Since we restrict ourselves to the consideration of the 
interaction regime where the dipole approximation ap- 
plies, we can further simplify the problem by reducing the 
dimensionality. We only treat the muon motion in one 
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dimension (ID) along the laser polarization axis. This 
is indicated, as in our case the ID numerics still is a 
non-trivial task because of the fine grid spacing required 
to resolve the nuclear extension and the non-standard 
laser parameters employed. As regards ordinary atoms, 
the latter would correspond to intense fields in the mid- 
infrared, giving rise to large muon momenta and pondero- 
motive energies. ID models are widely used in strong- 
field physics ^ as they retain the essential physical fea- 
tures of nonrelativistic laser-atom interaction for linearly 
polarized fields [ssj . 

The main shortcoming of the ID approach is the ne- 
glect of the muon's wave packet spreading in transver- 
sal direction. During one laser period T = 211/10, the 
spreading can be estimated as Ar AvT, where Av ~ 
(?ie_E/-^/2mp7p)^/^ is the velocity width at the exit of the 
potential barrier which the muon tunnels through [23| 
and Ip denotes the atomic ionization potential. In an ul- 
trastrong VUV field (/ = 10^3 W/cm^, huj = 60 eV) we 
obtain Ar ~ 10 pm for muonic hydrogen, which is by a 
factor p smaller than the spreading of an electron wave 
packet in an appropriately scaled infrared laser field [cf. 
Eq. (O]. The wave packet spreading can substantially re- 
duce the total harmonic yield. We stress, however, that 
the goal of the present study is to reveal relative differ- 
ences between physical observables which typically are 
less sensitive to model assumptions than absolute num- 
bers. 

In order to unveil the effects of the nuclear mass and 
size in the HHG spectra of muonic atoms we employ the 
following potentials in the ID version of Eq. ([3]); 



density within the nuclear radius R. The corresponding 
potential is 

v*W= l,^'- '"'^ (9) 

which explicitly takes the nuclear radius into account. 
We point out that in the limit R ^ 0, the binding energy 
of the lowest lying state of this potential becomes infinite 
and, thus, un phy sical [soj . Following a well-established 
procedure |39l44l| , we therefore start our calculation from 
the first excited state which has the correct binding en- 
ergy. By monitoring the projection onto the unphysical 
state during the time evolution, we take care that the 
occupation of this state always stays negligibly small. 



III. RESULTS 

The TDSE ([3]) for the relative motion has been solved 
numerically in one spatial dimension via the Crank- 
Nicolson time-propagation scheme. The laser field is al- 
ways chosen as a 5-cycle pulse of trapezoidal envelope 
having one cycle for linear turn-on and one for turn-off. 
The HHG spectrum is obtained from a Fourier transfor- 
mation of the dipole acceleration. 



A. Maximum cutoff energies 



1. Soft-core potential 

Numerical calculations in reduced dimensionality usu- 
ally employ a soft-core potential to describe the muon- 
nucleus interaction [ssj. In this manner, the Coulomb 
singularity of a pointlike nucleus at the origin is avoided. 
We also apply a standard soft-core potential in our com- 
putations, which after appropriate scaling reads 



Vsix) = - 



(8) 



It enables us to study the influence of the nuclear mass, 
which enters Eq. ([8]) via the reduced mass contained in 
the parameter p. The corresponding results are shown in 
Sec.IH.B. 



2. Hard-core potential 

In order to describe the effect of the finite nuclear ex- 
tension, the softcore potential ([8]) is not suitable. In- 
stead, we apply for this purpose the nuclear drop model 
and consider the nucleus as a sphere of uniform charge 



Since the conversion efficiency into high harmonics is 
rather low 10^^), it is generally desirable to maximize 
the radiative signal strength. In our situation, the opti- 
mization is of particular importance as the target density 
of muonic atoms is low. A sizeable HHG signal requires 
efficient ionization on the one hand, as well as efficient re- 
combination on the other hand. The former is guaranteed 
if the laser peak field strength lies just below the border 
of over-barrier ionization (OBI) where the Coulomb bar- 
rier is suppressed all the way to the bound energy level 
by the laser field [2l| . From Eq. ([3|) we obtain 
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with the fine-structure constant a « 1/137. Efficient 
recoUision is guaranteed if the magnetic drift along the 
laser propagation direction can be ignored, which limits 
the relativistic parameter to [13, El] 
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The condition (ITT]) also confirms the applicability of the 
dipole approximation in Eq. 
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The Eqs. pO|) and above define a maximum laser 
intensity /max ~ /'-'^i and a minimum laser frequency 
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(12) 



which are still in accordance with the conditions im- 
posed. At these laser parameters, the maximum har- 
monic cutoff energies are attained, while Eqs. ([T0|) and 
(|lip guarantee an efficient ionization-recollision process. 
Note, however, that smaller driving frequencies generally 
lead to reduced harmonic signal strengths because of the 
more pronounced and unavoidable quantum wave-packet 
spreading [43 - l46j (for an exeption to this rule, see 'i^]). 
For muonic hydrogen the lowest frequency according to 
Eq. (|12l) lies in the VUV range, fu^min ~ 27 eV, while the 
maximum field intensity is /°^^ « 1.6 x 10^'^ W/cm^. At 
these values, the harmonic spectrum extends to a maxi- 
mum energy of emax ~ 0.55 MeV. For light muonic atoms 
with nuclear charge number Z > 1, the achievable cutoff 
frequencies are even higher, reaching several Me Vs. A 
summary is given in TableUl 

For comparison we note that the highest harmonic cut- 
off energy which has been attained experimentally with 
ordinary (helium) atoms, amounts to « IkeV [48|; the 
corresponding harmonic order at the cutoff was emax/i^ ~ 
800. Higher cutoff energies are difficult to achieve due to 
the detrimental effects of dephasing '4^ and electron drift 
motion; various schemes have been proposed to overcome 
this obstacle (see [U H^l and references therein). 
Muonic atoms are advantageous in this respect since the 
large muon mass in principle allows for the generation of 
MeV harmonics in the dipole regime of interaction. 

As a result, muonic atoms are promising candidates for 
the generation of hard x-rays or even 7-rays which might 
be employed to trigger photo-nuclear reactions. 



z 




^min 


^max 


tmax 


1 


27 eV 


0.007 


0.085 


0.55 MeV 


2 


170 eV 


0.015 


0.12 


1.1 MeV 


4 


960 eV 


0.03 


0.17 


2.2 MeV 


10 


9.5 keV 


0.07 


0.27 


5.7 MeV 



TABLE I: Maximum HHG cutoff GIlGrglSS 6max achievable 
with hydrogenlike muonic atoms of nuclear charge number 
Z. The applied laser frequency ^min and intensity parame- 
ter ^max are chosen in accordance with Eqs. (|10|l - (|12| l to allow 
for an efficient ionization-recollision process, ^min denotes the 
minimum intensity parameter leading to tunneling ionization 

m. 



B. Nuclear mass effects 

In this section we consider the effect on the HHG pro- 
cess stemming from a variation of the nuclear mass, as- 
suming the nucleus as being point-like. We have solved 
the Schrodinger equation ([3]) using the softcore potential 



([8]). The latter depends on the nuclear mass via the re- 
duced mass entering the parameter p. It is meaningful 
to consider light isotopes where the reduced mass is 
significantly different from the bare muon mass . The 
nuclei chosen for the calculations in this and the follow- 
ing section are given in Tablellll In order to demonstrate 
the nuclear mass effect on the HHG process, isoptopes of 
large relative mass difference have been selected. 



Isotope 


Mass m„ (GeV/c^) 


Ref. 


Size R (fm) 


Ref. 


H 


0.9383 


[5^ 


0.875 


[52] 


D 


1.8756 


[5^ 


2.139 


[52] 


^He 


2.8084 


[5^ 


1.9448 


[53] 


^He 


3.7274 


[5^ 


1.6757 


[53] 


«Li 


5.6016 


[5^ 


2.517 


M 


''Li 


8.4069 


[5^ 


2.217 


M 


2«Ne 


18.493 


[55] 


3.0053 


[53] 


2^Ne 


21.277 


[55] 


2.9126 


[53] 



TABLE II: Nuclear masses and rms charge radii of various 
isotopes which have been used in the calculations. The half 
lives of ^Li and '^'^Ne amount to 178 ms and 37 s, respectively; 
the other nuclei are stable. 

First of all let us consider muonic hydrogen isotopes 
exposed to very intense VUV laser fields. In Fig.[lja) we 
show the harmonic spectra for muonic hydrogen (where 
the nucleus is a proton) and muonic deuterium. For 
muonic hydrogen the spectrum extends further including 
60 more harmonics as compared with that of deuterium. 
The difference of cutoff positions can be understood by 
inspection of the formula for the spectral cutoff energy, 
fmax = Ip + 'i.lTUp [SS]. In the present case the pondero- 
motive energy is given by [see Eq. (|3|)] 



1 

m,. 



1 



(13) 



0.90TOp, whereas m 



(D) 



and is, thus, the larger the smaller the reduced mass is. 
Consequently, in an intense laser field with Up 3> Ip, 
muonic hydrogen (H) gives rise to a larger cutoff energy 
than muonic deuterium (D). The relative difference is 
about 5% according to e^iSx/enSx ~ m^r^^ /m^^ 1.05. 
Note that TOr^-* « 0.90™^, whereas mi^'^ « 0.95mp 

The nuclear mass effect can also be understood more 
intuitively within the two-particle picture, instead of the 
relative motion. The right-hand side of Eq. (IT^ describes 
a ponderomotive energy that consists of two parts: one 
for the recolhding muon and one for the recolliding nu- 
cleus. The total pondermotive energy is thus a sum of 
the ponderomotive energies of the muon and the nucleus. 
Both particles are driven into opposite directions by the 
laser field and when they recolhde, their kinetic energies 
add up. In this picture the higher cutoff energy of the 
hydrogen atom arises from the larger ponderomotive en- 
ergy of the proton as compared to the heavier deuteron 
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In the opposite situation when Up <C Ip, the order 
of the spectral cutoff positions is reversed, as shown in 
Fig.[ljb) . Here we assumed a laser field with the same fre- 
quency as in Fig.[TJa) but a 100 times reduced intensity. 
The spectra in Fig.[IJb) exhibit a perturbative decay of 
the signal strength with increasing harmonic order, fol- 
lowed by a large peak around the 35th harmonic. The 
peak is due to a multiphoton resonance with the first 
excited atomic state. For muonic hydrogen, the transi- 
tion energy to this level is 2.00 keV in the softcore bind- 
ing potential © HI], whereas it amounts to 2.10 keV 
for muonic deuterium. When scaled to the laser pho- 
ton energy, the peaks arise at the corresponding, slightly 
different positions. 

For the parameters assumed in Fig.flja), the ratio 
fkd/Ip w 0.024 is relatively low. In fact, the same ratio is 
obtained when an ordinary hydrogen atom interacts with 
a mid-infrared laser field of frequency huj = 0.33 eV (i.e. 
wavelength A — 3.7 /im). Under such circumstances the 
harmonic signal strength is significantly reduced as com- 
pared to HHG in optical or near- infrared fields j43 - |46j . 
For this reason we provide in Fig.[ljc) the results for a 
higher driving frequency. As compared to the value em- 
ployed in Fig.ITJa), the laser frequency has been doubled 
leading to an increase of the harmonic signal strength 
by 1-2 orders of magnitude. Since the harmonic order 
at the cutoff position is proportional to uj^'^, the range 
of emitted harmonics is reduced. At the doubled laser 
frequency in Fig.[IJc), the different cutoff positions for 
muonic hydrogen vs. deuterium due to the nuclear mass 
effect are still clearly visible. However, at the tripled 
laser frequency (not shown) the cutoff positions practi- 
cally coincide because the cutoff energies from muonic 
hydrogen and deuterium differ by less than a photon en- 
ergy. We point out that the spectra can be distinguished 
nevertheless by the plateau height which remains sub- 
stantially higher for muonic hydrogen [see Figs.[lla) and 
He)]. 

In Fig. [5] the harmonic spectra for different isotopes 
of muonic helium in an ultra-intense XUV laser field are 
shown. We see that for ■^He the spectrum extends slightly 
further including 10 more harmonics in contrast with that 
of '*He. The difference in the cutoff positions cannot fully 
be explained by the reduced muon mass here. Formula 
which holds for hydrogen atoms {Z = 1), would 
predict a difference of only 5 harmonics for this case. In 
the general case {Z > 1), however, the ponderomotive 
energy of the relative motion reads 



The reason is that in the case Z > I, the center-of-mass 
does not stay at rest. Rather, the ponderomotive energy 



„2 Ri2 



1 



(14) 



This formula correctly predicts the difference in the cutoff 
positions in Fig. [2] The effective charge can thus have a 
measurable impact on the HHG response. 

It is interesting to observe that Eq. ([M]), in contrast to 
Eq. (jl3p , does not simply separate into a sum of the pon- 
deromotive energies of the muon and the nucleus; i.e., 
it is different from U'^ = [e^ E"^ j Au?){Z'^ j + \|m^,). 



(cm) 



U'p 

tion 



{Z - \)'^ I A^up- M is connected with its mo- 
The relation between the various ponderomotive 
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FIG. 1: (color online). HHG spectra calculated with the soft- 
core potential The black and grey (red) lines represent 
the spectra for muonic hydrogen and deuterium, respectively, 
(a) The laser parameters are / = 1.05 x 10^'^ W/cm^ and 
fej = 59 eV, corresponding to ^ ~ 0.04. (b) Same as (a) but 
for the reduced laser intensity I = 1.05 x 10^^ W/cm^. (c) 
Same as (a) but for the doubled laser frequency &j = 118 eV. 
Note that (a) has been reproduced from Ref. [17| in order to 
facilitate the comparison between the various spectra. 
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Harmonic Order 



FIG. 2: (color online). HHG spectra calculated with the soft- 
core potential (|8]). The black and grey (red) lines represent 
the spectrum for muonic "^He and ''He, respectively. The laser 
parameters are I = 8 x lO^"* W/cm^ and tiw = 347 eV. The 
inset shows an enlargement of the cutoff region. 

energies is 

t/; = c/w + f/(™). (15) 

Only in the case of hydrogen isotopes {Z ~ 1) the center- 
of-mass coordinate remains at rest since the total charge 
is zero, so that Up'^™^ = and Up^^ fully accomodates the 
single-particle ponderomotive energies. 

C. Nuclear size effects 

In this section we investigate the influence of the nu- 
clear extension on the radiation spectra by making use 
of the hardcore potential (jH]). For the radial parameter 
R in Eq. ([9]) , the respective rms radius from TableHl] is 
employed. We again consider low-Z isotopes where the 
largest relative size differences are found. Two isotopes 
with different radius also differ in mass. We wish to sep- 
arate, though, the impact of the nuclear size from the 
nuclear mass effect which was discussed in the previous 
section. To this end, we compare in the following HHG 
spectra from different isotopes where the nuclear mass 
effect has been removed by a suitable adjustment of the 
laser frequencies and intensities. In accordance with the 
scaling relations (I?]), this was achieved by applying the 
scaled parameters lu cx rrir and E oc m^/qe (at a given 
value of Z). In this manner, the laser-driven muonic iso- 
topes become equivalent to the same ordinary hydrogen 
atom, with the only difference being the size of the bind- 
ing nucleus. In particular, the harmonic cutoff positions 
are forced to coincide this way. 

We start with the calculations for muonic hydrogen 
and deuterium shown by Fig. [31 As mentioned above, 
in order to avoid residual signatures from the nuclear 
mass effect we apply in Fig.[3][a) the laser parameters 
/(H) = 1.05 X 1023w/cm2 and = 59 eV to muonic 

hydrogen, whereas muonic deuterium is subject to the 




Harmonic Order 




Harmonic Order 



FIG. 3: (color online). HHG spectra calculated with the 
hard-core potential ((9|. (a) The black line represents the 
spectrum for muonic hydrogen at the laser parameters 7'^' = 
1.05 X 10^3 W/cm^ and Hlj^^^ = 59 eV. The grey (red) line 
represents the spectrum for muonic deuterium at the appro- 
priate scaled laser parameters l'-^'' = 1.30 x lO^^W/cm^ 
and = 62 eV in order to compensate for the nuclear 

mass effect. The inset shows an enlargement of the cutoff 
region on a linear scale (reproduced from [l7j to allow for a 
direct comparison), (b) Same as (a) but for the laser pa- 
rameters = 1.05 X lO^^W/cm^ huu'-^'' = 177 eV and 
1^°^ = 1.30 X lO^^W/cm^, fkj^°'> = 186 eV, whereas (c) 
refers to = 3.8 x 10^^ W/cm^ Hlj'-^'' = 59eV and 
J(°' = 4.68 X 10^^ W/cm^ huj'^^^ = 62 eV. 
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parameters = 1.30x10^^ W/cm^ and fia;(°) = 62 eV. 
Regarding the overall shape of both spectra we observe 
a dip at low harmonics around the order n ^ Ip/Huj, fol- 
lowed by a rising plateau region. These common features 
are due to the hardcore potential used, in contrast to 
the softcore potential ([8|). In fact, a siinilar behaviour of 
the harmonic response was found in [57] and attributed 
to so-called non-tunneling harmonics in very steep poten- 
tials. Our main focus, however, lies on the relative differ- 
ence between the two spectra. The harmonic signal from 
muonic hydrogen is larger (by about 50 % in the cutoff 
region) than that from muonic deuterium. The reason is 
that in the case of muonic hydrogen the nuclear radius is 
smaller which generates a steeper potential near the ori- 
gin [see Eq. ([9|)]. This leads to a larger potential gradient 
in this region which accelerates the atomic dipole accord- 
ing to Ehrenfest's theorem j58| . The muon in hydrogen 
is thus more strongly accelerated, leading to enhanced 
harmonic emission. 

When the laser frequency is enhanced three times as 
compared to Fig.[3l[a) while the laser intensity is kept 
constant, the relative difference between the HHG spec- 
tra from muonic hydrogen and deuterium stays approxi- 
mately constant. As shown in Fig. [3Kb), the signal from 
muonic hydrogen remains larger by about 50 % in the 
cutoff region. When instead the frequency of Fig.jSja) is 
kept but the field intensity is decreased down to 3.8 x 10^^ 
W/cm^ (4.68 X 10^^ W/cm^) for muonic hydrogen (deu- 
terium), then the harmonic signals only differ by about 
10% in the cutoff region, as shown in Fig.[3fc). 

In the case of muonic helium isotopes, the harmonic 
signal from "^He is expected to be larger than that from 
•^He. For *He, though containing one more neutron, is 
a doubly magic nucleus of very compact size (see Ta- 
ble lll)) . This expectation is confirmed by Fig.UJ^a), where 
a relative difference of about 10 % in the cutoff region 
is observed. The fact that the difference is reduced in 
comparison with the muonic hydrogen isotopes can be 
attributed to the smaller relative difference of the nuclear 
radii: i?('"<^)/i?(''H<=) w 1.16, whereas R^^^R^^^ « 2.44. 
On the other hand, however, the muon comes closer to 
the binding nucleus when the charge number Z increases 
which should enhance the sensitivity to the nuclear size. 
This circumstance becomes important when we move on 
to muonic lithium (Z = 3). Here, the diflterence between 
the harmonic signals from muonic ^Li versus muonic ^Li 
amounts to about 20 % in the cutoff region [see Fig.|4l[b)] 
which is larger than for the helium isotopes although the 
ratio of the nuclear radii R^ ^'^i?^ '^'^ w 1.14 is simi- 
lar here. In order to facilitate the comparison between 
helium and lithium, the laser parameters in Fig. 3] were 
chosen to generate a uniform cutoff position. Finally, a 
comparison of the radiative responses from muonic ^"^Ne 
and ^"^Ne at the corresponding field parameters reveals 
almost identical HHG spectra (not shown). The relative 
enhancement of the HHG signal from the smaller isotope 
^^Ne as compared with ^'^Ne is of the order of 1%. Here, 
the ratio of the nuclear radii, i?("'N<=V-R^''^°-' ~ l-03> 
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FIG. 4: (color online) . HHG spectra calculated with the hard- 
core potential ((9)|. (a) The black line shows the spectrum for 

muonic ''He at the laser parameters /' = 8.7 x 10^'' W/cm^ 
and fioi' — 255 eV. The grey (red) line represents the 
spectrum for muonic ^He at the accordingly scaled values 
jr(3Hc) ^ g g ^ io24w/cm2 and tnu^"'^'''' = 253 eV. (b) Same 
as (a) but for muonic ^Li at /'"^'^ = 1.06 x 102''W/cm^ 
fej'^"^'' = 583 eV [black line] and muonic '^Li at j'"'^'' = 
1.01 X lO^** W/cm^, fiajf^Li) = 580 eV [grey (red) line]. 



is close to unity. We note moreover that in the cases 
of lithium and neon practically no mass effect needs to 
be compensated since the corresponding reduced muon 
masses almost coincide. 

So far we have restricted our investigation of the nu- 
clear size effects to isotopes which exist in nature. It 
is instructive to extend the consideration by varying the 
nuclear radius artificially in order to gain a more com- 
plete picture of the physical mechanisms involved. Con- 
sidering hydrogenic nuclei with charge number Z = 1 
and radius R increasing from 1 fm to 5 fm, we obtain 
a monotonously decreasing HHG signal, in accordance 
with the explanation given above regarding Fig. [3] For 
neon-like nuclei with Z = 10, however, a minimum har- 
monic emission strength arises around i? ~ 3 fm. When 
the nuclear extension is artificially enhanced further to 
R = 5fm, the HHG signal is growing again (see Fig. [5]). 
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R{fm) 

FIG. 5: Dependence of the harmonic signal strength at the 
cutofT on the radius of the neon-hke nucleus to which the muon 
is bound. The signal strength is expressed as the normalized 
difference A = [S{R)-S{Ro)]/S{Ro), where SiR) denotes the 
near-cutoff signal strength for nuclear radius R and i?o = 3 fm 
is the reference value. We point out that the nuclear size 
variation shown here largely extends beyond the Ne isotopes 
which exist in nature where always _R ~ 3 fm 53] . 



This different behavior can be attributed to the atomic 
ionization potential Ip which significantly decreases in 
the neon-like case when R in Eq. ([9]) is increased. Note 
that the atomic Bohr radius of hydrogenlike muonic neon 
amounts to 25 fm only, rendering this system more sen- 
sitive to strong deviations from a point nucleus. Hence, 
the probability for tunneling ionization in the laser field, 
being the first step of HHG, is enhanced which leads to 
amplified harmonic emission for R > 3fm. 

When considering laser-driven recollisions, the de- 
Broglie wavelength of the returning quantum wave packet 
can be of importance as well. In fact, for the case of 
ordinary molecules it has been shown that the electron 
wavelength can become as small as the internuclear dis- 
tance within the molecule, causing characteristic diffrac- 
tion patterns ^59^. In the present case of muonic atoms, a 
similar effect could in principle arise when the de-Broglie 
wavelength of the recolliding muon compares with the nu- 
clear size. However, the wavelength of a muon with a ki- 
netic energy of a few MeV (see Table|l]) amounts to about 
50 fm which exceeds any nuclear radius substantially and 
thus prevents diffractive muon-nucleus scattering. 

Concluding this section, we have shown that the 
plateau height of HHG spectra from muonic atoms is 
sensitive to the finite nuclear size. Smaller nuclei within 
the range of existing isotopes lead to enhanced harmonic 
emission. We note, however, that - contrary to the nu- 
clear mass effect in Sec. III.B - the influence of the nu- 
clear size might be overestimated by our ID hardcore- 
potential approach as the muon meets the nucleus more 
often than in the real 3D case. A calculation in higher 
dimensionality could provide more accurate quantitative 
predictions on the nuclear size effect whose physical ori- 
gin and basic features have been presented here. 



D. Comparison with electronic systems 

Finite nuclear size effects - in the absence of any ex- 
ternal laser field - have been revealed in high-precision 
spectroscopy of electron transitions in (ordinary) highly 
charged ions (see, e.g., [gO, HH for recent experiments). 
When such ionic systems are exposed to a superintense 
laser field , nuclear signatures may be present in their 
high-harmonic response as well. It is of interest to com- 
pare the expected effects with those found for laser-driven 
muonic atoms in Sec. III.C. To this end we perform a sim- 
ple analysis which is based on nonrelativistic Schrodinger 
theory; the relativistic electron motion in highly charged 
ions is ignored in this rather qualitative discussion. 

We assume a hydrogenlike system of nuclear charge 
number Z and employ the mass scaling parameter p 
from Eq. (5), with p « 1/200 for a muonic atom 
and p = 1 for an electronic ion. The K-shell Bohr 
radius, binding energy, and Coulombic field strength 
amount to aK{Z,p) — aop/Z, Ip{Z,p) = eaZ'^/p, and 
Ek{Z,p) = EqZ^/p^, respectively, where ao, eo and i?o 
denote the corresponding quantities for ordinary hydro- 
gen. The nuclear radius can be approximated roughly 
as R{Z) « 1.2(2Z)^/'^ fm and has a typical relative vari- 
ation among different isotopes of a few percent (except 
for hydrogen vs. deuterium). Similar finite nuclear size 
effects in the HHG spectra can be expected when the ra- 
tio R{Z)/aK{Z, p) oc Z'^l'^ I p has a similar value for two 
atomic systems that are compared. This is the case, e.g., 
for electronic U^^+ (where Z = 92, p = 1) and muonic 
He+ (where Z = 2,p w 1/200). 

The above relations imply, however, that the binding 
energy and electric field strength in the electronic ion 
are substantially larger than in the muonic atom when 
both have the same ratio of Z'^l'^ j p — const. As a con- 
sequence, the laser frequency and intensity that must be 
applied to the electronic highly charged ion in order to 
reveal finite nuclear size effects in the harmonic response, 
need to be larger than in the muonic atom case. Against 
this background, muonic atoms appear as more favorable 
systems than ordinary heavy ions to study the influence 
of the nuclear size on the HHG process. 



IV. CONCLUSION AND OUTLOOK 

Motivated by the sustained progress in the develop- 
ment of powerful laser sources, we have studied the har- 
monic radiation which is emitted by muonic atoms ex- 
posed to high-intensity, high-frequency laser flelds. It 
was shown that maximum harmonic cutoff energies in the 
MeV domain can be achieved, rendering this species of 
exotic atoms promising candidates for the generation of 
(weak) ultrashort coherent 7-ray pulses which might be 
employed to trigger photo-nuclear reactions. Our results 
demonstrate moreover that strongly laser-driven muonic 
atoms can, in principle, be utilized to dynamically gain 
structure information on nuclear ground states via their 
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high-harmonic response. (1) On the one hand, the har- 
monic cutoff position extends to larger values for isotopes 
of smaller mass. For the hydrogen isotopes this effect is 
fully explicable in terms of the reduced mass, whereas for 
atomic numbers Z > 1 an additional contribution stems 
from an effective muon charge which affects the relative 
motion generating the harmonics. (2) On the other hand, 
the harmonic signal strength additionally depends on the 
nuclear size, being enhanced for more compact isotopes. 
Corresponding nuclear size effects in the high-harmonic 
emission from ordinary highly-charged ions are expected 
to be less pronounced. 

Furthermore, we point out that the interaction of a 
muonic atom with ultrastrong laser fields may lead to 
excitation of the nucleus. Nonresonant nuclear Coulomb 
excitation has recently been studied when a bound 
muonic wave packet is driven into coherent oscillations by 
an external laser field [gS]; the resulting nuclear excita- 
tion probabilities were found to be small, though. When 
the laser field is sufficiently strong to ionize the muon as 
in the HHG scenario, however, the kinetic energy gain in 



the continuum up to the MeV range allows for nuclear ex- 
citation upon the muon-nucleus recollision. Correspond- 
ing studies of laser-driven electron-impact excitation of 
the nucleus have been carried out in ordinary atoms and 
ions [13] ■ It is even conceivable to conduct pump-probe 
experiments on excited nuclear levels: the periodically 
driven muon can first excite the nucleus and then probe 
the excited state and its deexcitation mechanism during 
a subsequent encounter. 
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